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Single crystal superalloys have become increasingly popular for turbine blade and 
vane applications due to their high strength, and creep and fatigue resistance at elevated 
temperatures. The crystallographic orientation of a single crystal material greatly affects 
its material properties, including elastic modulus, shear modulus, and ductility. These 
directional properties, along with the type of loading and temperature, dictate an 
anisotropic response in the yield strength, creep resistance, creep rupture ductility, 
fatigue resistance, etc.. A significant amount of research has been conducted to 
determine the material properties in the <00 1> orientation, yet the material properties 
deviating from the <00 1> orientation have not been assessed for all cases. Based on the 
desired application and design criteria, a crystal orientation is selected to yield the 
maximum properties. Currently, single crystal manufacturing is able to control the 



primary crystallographic orientation within 15° of the target orientation, which is an 
acceptable deviation to meet both performance and cost guidelines; the secondary 
orientation is rarely specified. 

A common experiment is the standard load-controlled tensile test, in which 
specimens with different orientations can be loaded to observe the material response. 

The deformation behavior of single-crystal materials under tension and compression is 
known to be a function of not only material orientation, but also of varying micro- 
deformation (i.e. dislocation) mechanisms. The underlying dislocation motion causes 
deformation via slip, and affects the activation of specific slip systems based on load and 
orientation. The slip can be analyzed by observing the visible traces left on the surface of 
the specimen from the slip activity within the single crystal material. The goal of this 
thesis was to predict the slip systems activated in three-dimensional stress fields of a 
notched tensile specimen, as a function of crystal orientation, using finite element 
analysis without addressing microstructural deformation mechanisms that govern their 
activation. Out of three orientations tested, the specimen with a [1 10] load orientation 
and a [001] growth direction had the lowest maximum resolved shear stress; this 
specimen orientation appears to be the best design candidate for a tensile application. 


CHAPTER 1 
INTRODUCTION 


Project Background and Goals 


Single Crystal Superalloys 

Nickel-base single crystal superalloys have become increasingly popular, 
particularly in turbine blade applications, because of their exceptional thermo mechanical 
fatigue properties at high temperature. Unlike the more commonly used isotropic alloys, 
these superalloys are orthotropic and have highly directional material properties because 
they are grown as single crystals (i.e. as one grain). The most common primary growth 
direction for the nickel-base superalloys is the <00 1 > direction; <00 1 > is not only the 
most easily grown, but is also the direction with the most desirable combined strength 
properties. This is advantageous because many parts are cast, or grown, rather than 
manufactured from a larger single crystal sample with a specific orientation. Currently 
we are most interested in one particular alloy, which will be referred to here as “Material 
A,” because of its use in advanced aircraft engines. (Note: Currently we are unable to 
disclose the specific alloy due to security restrictions.) These superalloys play an 
important role in commercial and military propulsion systems as well. 

Single crystals are manufactured by selecting grains with a desired orientation 
and then “growing” the part in one of two ways. The first method uses a helical mold, 
commonly called a “pigtail,” placed between a chill plate and the part casting. As a grain 


moves through the mold, the helix enables selection of grains with a <00 1 > orientation, 
which then grow through the coils and up through the mold itself. Although some 



Figure 1-1 Helical mold cast turbine blade. 

Source: Deluca, 2001. 

undesirable grains may initially form at the starter plate, after building up through the 
coils only the <00 1 > grain orientation remains (Figure 1-1). The coils are later removed 
to leave the desired part geometry. The helical mold is frequently used to manufacture 
turbine blades, but it can only be used to cast parts with a primary orientation along one 
of the <100> directions due to the helix restrictions combined with the material’s 
orthogonal grain growth (Davis, 1997). The second method uses a grain starter, or seed, 
and is consequently termed “seeding.” Like the helical mold, this method can also be 
used to grow <100> primary orientations, but it can also produce <01 1> and <1 1 1> 
crystals as well. In fact, a properly chosen seed grain can produce any desired primary 


and secondary orientation combination. The seed is placed on a chill plate at the base of 
a casting filled with the molten material. The temperature is tightly controlled to prevent 
the seed from melting, and the alloy in the mold conforms to the seed alignment as the 
material cools. 

Currently, the primary direction (i.e. growth direction) can be controlled 
relatively well, but still has an allowable deviation of 1 5° in any direction. Although it is 
possible to reduce this deviation by simply discarding those parts with larger variations, 
the cost in wasted materials would be too great and would result in a significant 
percentage of scrapped parts. Nonetheless, the material properties vary greatly as the 
primary direction deviates from the <00 1 > orientation, and it is important to quantify the 
effect of other 


Airfoil 

Stacking 

Line 



Figure 1-2 Primary and secondary turbine blade crystallographic orientations. 

Source: Moroso, 1999. 



orientations on material properties. The secondary orientation can also be controlled 
during casting (Figure 1-2), however most manufacturers choose to ignore this control to 
achieve greater productivity. The material can, however, be examined after casting to 
determine secondary orientation, and if a particular secondary orientation has beneficial 
material properties for the specific application the part can then be further manufactured 
by cutting from the cast piece. When it becomes possible to more tightly control the 
manufacturing process, both the primary and secondary orientations with the best 
material properties can be combined and applied. Although <00 1 > is well known to be 
the strongest orientation in terms of yield strength (Davis, 1997), the underlying 
dislocation and slip deformation mechanisms, with respect to the primary orientation, are 
still under study. 

Sli p Deformation 

Nickel-base superalloys have a basic face centered cubic (FCC) structure in the matrix 
phase. These superalloys undergo the usual method of plastic deformation by slip, where 
separate parts of the crystal structure slide over one another along definite 
crystallographic, or slip, planes in specific directions (Dieter, 1986). When a load is 
applied to an FCC single crystal specimen, the first slip systems to be activated are those 
termed the “easy glide” systems (Table 1-1), or the primary octahedral slip systems. 
Typically, these are the planes with the greatest atomic densities. As deformation and 
slip activity continue, the material begins to exhibit physical evidence of this slip at the 
surface (Figure 1-3). Slip occurs so that the energy of the high shear stress within these 
12 primary slip systems (knows as the “resolved shear stress” or RSS) is alleviated; the 



Table 1-1 Slip systems in an FCC crystal. 


Slip Number 

Slip Plane 

Slip Direction 


Octahedral Slip a/2{1 1 1}<1 10> 

1 

(111) 

[10-1] 

2 

(in) 

[0-11] 

3 

(HD 

[1-10] 

4 

(-ii-i) 

(10-1] 

5 

(-11-D 

[110] 

6 

(-11-1) 

[011] 

7 

(1-1-D 

[110] 

8 

(1-1-D 

[0-11] 

9 

(1-1-1) 

[101] 

10 

(-1-11) 

[011] 

11 

(-1-11) 

[101] 

12 

(-1-id _ 

U-io] 


Octahedral Slip a/2{1 1 1}<1 12> 

13 

(in) 

[-12-1] 

14 

(in) 

[2-1-1] 

15 

(in) 

[-1-12] 

16 

(-11-1) 

[121] 

17 

(-ii-i) 

[1-1-2] 

18 

(-ii-i) 

[-2-11] 

19 

(i-i-i) 

[-11-2] 

20 

(i-i-i) 

[211] 

21 

(1-1-1) 

[-1-21] 

22 

(-i-ii) 

[-21-1] 

23 

(-i-ii) 

[1-2-1] 

24 

(-i-ii) 

Q12] 


Cubic Slip a/2{100)<1 10> 

25 

(100) 

[011] 

26 

(100) 

[01-1] 

27 

(010) 

[101] 

28 

(010) 

[10-1] 

29 

(001) 

[110] 

30 

(001) 

[-110] 


Source: Stouffer and Dame, 1996. 



Figure 1-3 Slip lines on the surface of a [100] loaded tensile specimen. 

Source: Forero and Ebrahimi, 2002. 

slip, or dislocation motion, results in a more energetically stable system. Therefore, 
materials typically fail along one of the primary {111} octahedral planes (Figure 1-4). 
As the applied stress continues to rise, or with changes such as loading type/direction, 
time or temperature, the activation of 12 secondary slip systems may occur. At low 



Figure 1-4 Turbine blade failure on a {1 1 1} octahedral plane. This particular turbine blade 
failed because the attached machine’s frequency was equivalent to the turbine blade’s natural 

frequency (upper right comer). 

Source: Deluca, 200 1 . 


modulus orientations, or at extreme temperature/load conditions, the initiation of six 
cubic slip systems may also occur (Table 1-1). These secondary and cubic systems are 
not typically activated initially, because slip along the primary systems requires less 
energy to alleviate the high RSS. In addition to load and temperature conditions, the shift 
in slip from the primary to secondary or cubic slip systems is based on specific 
microstructural behaviors, such as the pinning or locking of dislocations, and will be 
discussed later only at a cursory level. 

Slip will occur when the RSS exceeds the yield strength of the material. 

However, nickel-base superalloys have been shown to exhibit different yield strengths in 
tension versus compression, a behavior known as tension/compression asymmetry. The 
<00 1 > orientation exhibits its highest yield strength in tension, while its compressive 
yield strength is lower. Conversely, a sample loaded in the <1 10> direction is stronger in 
compression rather than tension, while the <1 1 l>-loaded material has virtually no 
tension/compression asymmetry. Extended dislocations are dislocations whose parts 
separate to reach a lower energy state; these dislocations, as on octahedral slip planes in 
superalloys, must recombine into a single dislocation before cross-slip can occur. Cross- 
slip is the process where a dislocation moves from one plane to another, again to reduce 
the energy of the system. As Lall, et al. (1979) proposed, an applied load will either aid 
or hinder recombination of the dislocation parts. Since the material orientation affects 
the magnitude of the resolved shear stresses, it will also affect this process of 
recombination. For example, if a compressive load is applied and two partial 
dislocations are at an interface where cross-slip can occur, the applied load can overcome 
the force separating the two partials, and combine them to form a single dislocation. The 


unified dislocation can then move by cross-slip. If, however, a tensile load is applied to 
the same extended dislocation, the load will aid the force separating the partials and make 
cross-slip more difficult. 

Asymmetrical load behavior becomes extremely important when designing a part 
that must accept both tensile and compressive loads, for example parts under cyclical 
loading. However, the priority given to compressive stress analysis is a distant second 
behind tensile stress analysis with respect to fracture research, because compressive 
stresses are generally beneficial. Rather than causing cracks to initiate or propagate, 
compressive stresses usually either have no effect or may even arrest cracks that 
developed under tensile stresses. Therefore our focus on tensile testing is practical and 
relevant. 

Nickel-base superalloys exhibit another abnormal yield characteristic called 
“anomalous yield behavior.” Normal yield behavior shows decreasing yield strength 
with increasing temperature. As the temperature rises, diffusion enables cross-slip to 
occur more easily and the dislocations encounter fewer energy-based restrictions; the 
overall force required to produce deformation is reduced. The nickel-base superalloy 
matrix and precipitate show the opposite trend: as temperature increases their yield 
strength increases, up to a critical temperature. This behavior is a result of dislocation 
mechanisms and, like tension/compression asymmetry, is also affected by orientation. 

As the temperature increases and enables cross-slip to cubic planes, the dislocations 
dissociate to lower their energy and may form “locks,” (Kear-Wilsdorf locks, among 
others (Zhu, et al., 1998)). The term “lock” is used because the dislocation separates in 
such a way, on two separate planes, that no further motion can take place. Other 


dislocations then begin to build up at these locks, which act as barriers to further slip and 
raise the yield stress by impeding motion. Anomalous yield behavior is an important 
consideration for fracture and fatigue analysis, particularly at high temperatures, in order 
to understand mechanisms that increase strength and their limitations. For this behavior 
to occur, the RSS on the cubic planes will, presumably, be higher than the octahedral 
planes; if the cubic planes are taking more of the stress, and are able to contain it with the 
locks, slip will only be achieved when a sufficient stress is resolved onto the octahedral 
planes. If, however, the stress on the octahedral planes is higher it should result in the 
traditional yield behavior and “easy glide” slip regardless of the dislocations on the cubic 
planes. 

Our focus will be on the variation of the 12 primary resolved shear stresses in a 
notched specimen, and to determine the activity of the specific slip systems with respect 
to radial and angular distances from the notch tip for specific orientations (Figure 1-5). 


90 



p = notch radius 


Figure 1-5 Radial arcs for numerical stress field calculations. 



The slip systems will be examined both near the notch tip and at far-fleld to observe 
changes associated with the high stress gradients prevalent in close proximity to the 
notch tip. The variation of the RSS with respect to theta will show which systems are 
most active at various locations around the notch. What we expect to see, along a line of 
constant radius, is shifting maximum RSS values and slip systems, indicating a shift in 



Figure 1-6 Slip sectors occurring under plastic deformation. 

Source: Crone and Shield, 2001. 

the state of stress. These different slip systems can be clearly seen as “sectors” 
surrounding the notch tip in actual material tensile testing (Figure 1-6). The goal of this 
paper is to predict slip system activity, indicated by the maximum RSS, as a function of 
the radial and angular distances from the notch tip and the resulting slip sectors. 

Slip deformation is generally considered a stress-controlled process; the critical 
resolved shear stress (CRSS), or the stress at which the material is predicted to slip on 
any particular slip plane, is the controlling value. The CRSS is the shear stress that 


equals the material’s yield strength, and is a function of the applied load and direction, 
specimen geometry, and crystal structure. Therefore the RSS is a function of geometry, 
and the CRSS is a function of geometry and yield stress; neither are directly related to the 
material’s anisotropy. 

Material Characteristics 

The microstructure of nickel-base superalloys consists of a y-matrix and a fine 
dispersion of hard, y’-precipitates. The matrix, mainly nickel, is heavily alloyed with 
other elements that vary with the given alloy, including cobalt, chromium, tungsten and 
tantalum; the precipitate is the intermetallic compound N^Al. These superalloys have 
evolved in three “generations,” thus far. The most advanced, or third-generation, 
superalloys include Rene N4, CMSX-4, and others, and have a high volume fraction of 
y’, around 60% (Figure 1-7). (Note: The prefixes denote the company that manufactures 



Figure 1-7 Microstructure of Material A. The y’ precipitate (approximately 0.5gm in length) 
forms in the y matrix and comprises nearly 60% of the material. 

Source: Deluca and Annis, 1995. 



the material; e.g. CMSX = Cannon-Muskegon, where SX represents single crystal. The 
numbers are somewhat arbitrary depending on the particular manufacturer.) Although 
the majority of deformation occurs in the softer matrix (Svoboda and Lukas, 1998), the 
precipitate has a pronounced effect on the superalloy’s overall performance, especially at 
such high volume fractions. 

Temperature Effects 

One of the primary reasons for the popularity of Material A and other nickel-base 
superalloys is their high temperature capabilities. As noted previously, dislocation 
mechanisms can change with temperature, and a significant amount of research has been 
conducted on their behavior. Most notched tensile tests conducted to study slip sectors 
have been performed at standard room temperature; room temperature is in the low 
temperature regime, well below the transition temperature for superalloys. As Stouffer 
and Dame (1996) report, octahedral plane slip deformation usually dominates, and 
theoretically exclusively controls, low temperature deformation. As the temperature rises 
the secondary planes may be activated, and beyond T > 600°C thermally activated cube 
slip begins and acts along with octahedral slip. As the temperature further increases 
beyond approximately 850°C, climb and cross-slip dominate deformation and the 
material’s strength dependence on orientation is greatly reduced. Miner, et al. (1986) 
also found that near <1 1 1> only cubic deformation was evident, while near <1 10> and 
<100> only octahedral slip was present. Our study will focus on <11 0> and <100> 
orientations, supporting our focus on octahedral slip. 



Test Methods 


Several analytical and numerical solutions have been developed for the elastic 
response of isotropic notch specimens in tension, particularly in the field of linear elastic 
fracture mechanics. However, the isotropic analytical models have been very difficult to 
develop into three-dimensional anisotropic models; the current solutions rely on many 
simplifications that give inaccurate results when compared to experimental data. 

However, the capabilities of three-dimensional finite element analysis (FEA) can account 
for the current limitations in the elastic models and enable a solution that should correlate 
well to actual experimental results. Both the numerical and experimental specimens may 
include notches, which act as very simplified cracks to model fracture behavior. The 
effect of anisotropy is of great interest because many materials, although designed for 
primary strength in one direction, must withstand multi-axial loading. 

Analytical Approach 

Analytical solutions have the ideal trait of providing an exact solution to a 
problem. However, very complex problems often do not have exact analytical solutions; 
many approaches represent a combination of theoretical and empirical solutions, or close 
approximations. For the notched tensile specimen, an analytical mechanics solution is 
not available (much less one that can account for material variations) that correctly 
predicts slip sectors or slip activity about the notch tip. The current analytical solutions 
incorporate many assumptions, including plane strain, which will be discussed in depth in 


Chapter 2. 



Numerical Approach 


Finite element analysis is able to account for gross isotropic material properties, 
such as modulus of elasticity and Poisson’s ratio, as well as their directional counterparts 
in anisotropic materials. However, FEA typically does not account for the 
microstructural properties that dictate yield strength, such as dislocation mechanisms, nor 
can it predict other microstructural behavior. Small-scale atomistic simulations are 
capable of predicting dislocation generation, interaction, etc., however it would be much 
too costly to analyze even a very small specimen with actual dimensions on the atomic 
level. Limiting the size for a reasonable atomistic simulation would distort the model to 
the degree that the validity of the results would be highly questionable. Therefore, FEA 
is the only feasible type of computer simulation currently available to model notched 
single crystal specimens. Because it neglects microstructural behavior, FEA can also 
determine what influence the specimen’s geometry and anisotropy have alone on material 
property behavior, without the additional consideration of atomic interactions. 

Experimental Approach 

Tensile testing has long been used to measure material properties such as stress- 
strain behavior, yield strength, etc. and is common in many materials testing laboratories. 
A widely used experimental test sample is the notched specimen, loaded to observe the 
effect on the material’s overall displacement, stress and strain fields, and of particular 
interest here, slip line deformation. The notch introduces a triaxial state of stress in the 
proximity of the notch, thus providing an environment to study slip system formation in 
three-dimensional stress fields. Future tensile tests will be carried out for the single 
crystal superalloy Material A to observe the effect of the load orientation on the active 



slip regions, or “sectors,” about the notch to compare with FEA results. Using a constant 
load several tests specimens with different crystallographic orientations will be tested to 
failure to observe the active slip planes. Although this is clearly a different procedure 
than the purely elastic response measured by FEA, the magnitude of the applied elastic 
stress is an indication of which planes will first allow plastic deformation. If the most 
highly-stressed planes in the elastic analysis do not correlate to the slip lines observed 
from the tensile specimens, it will support the influence of other dislocation mechanisms 


on fracture. 



CHAPTER 2 
LITERATURE REVIEW 


An analysis of the stress-field at a single crystal notch involves two main 
components: incorporating the correct component and RSS calculations, and then 
defining the resulting slip zones around the notch. For a solid (i.e. no notch) specimen, 
the stresses can be precisely calculated using an analytical approach. Obviously there are 
not sectors defined for this solid material. However, the more complex, and relatively 
unexplored, study of the stresses and sectors around a notch tip is of recent origin; despite 
the analytical and experimental investigations by researchers (Rice, 1987; Shield, 1996; 
Schulson and Xu, 1997; Crone and Shield, 2001) an accurate single crystal model is far 
from complete. 


Slip Activation and Deformation 

Slip in single crystals should theoretically occur on octahedral planes (Chapter 1). In an 
isotropic material the twelve primary slip systems (or fewer depending on orientation), 
should be activated simultaneously based on equal Schmid factors (Figure 2-1). The 
Schmid factor, m, is a function of the load orientation, the slip plane orientation and the 
slip direction: 

m = cos X cos c|> (2-1) 


^rss = m ^ 


( 2 - 2 ) 



Where o is the applied load, T rss is the RSS component in a given slip plane and direction, 
X is the angle between the direction of the applied load and the shear direction, and (p is 
the angle between the applied load and the normal to the slip plane (Figure 2-2). 


Primary; t 1 , t 2 , t 3 Primary- r 4 , t 5 , 



Primary: x r . t* x * 


Primary: t 11 , t ^ 2 



Figure 2-1 Primary resolved shear stress planes and directions. 
Source: Modified from Stouffer and Dame, 1996. 




Figure 2-3 Persistent slip bands in an FCC specimen. Slip lines, formed by octahedral slip 
system activity, are clearly visible as the specimen undergoes plastic deformation. 
Source: Deluca and Annis, 1995. 


Following Dieter (1986) the cosine of the angle between two directions [hi ki \\] and 
[h 2 k 2 b] can be found using the direction indices: 


h ] h2 + k j • k2 + lj b 


(2-3) 


cos0 = 

f ~ 2 ~2 ~2 f ~2 ~2 ~2 

yj h 1 + kl + h yj^2 + ^2 + h 

Since the CRSS is reached when the RSS is equal to the yield stress of the material, the 
slip systems with the highest Schmid factors will reach the CRSS first: This is the so- 
called u Schmid Law” (Figure 2-3). However, compliance to Schmid Law has been 
proven only with isotropic materials, and its correlation to single crystals is not yet 
known. Therefore, another method must be used to predict RSS values and slip 
activation for these anisotropic materials. 

Anisotropy of Elasticity 


Elasticity 

Elasticity is defined by specific elastic constants that relate to atomic strength and 
spacing (Dieter, 1986). Elasticity, therefore, varies with orientation. In any given 
direction the spacing between atoms in an FCC unit cube is different. For example, in 
Figure 2-4 the distance from A to B is ao/2, while the distance from A to C is a 0 V2/2. 



Figure 2-4 FCC crystal structure. 



Another way of expressing the relative spacing between atoms is by the atomic density, 
or number of atoms per unit area (Table 2-1). The term a 0 refers to the unit atomic 

spacing, therefore a 0 2 is the area of the {100} planes, etc. The greatest atomic density 

correlates to the least distance between atoms, and the most likely planes for slip. Planes 
such as the {111} planes are called “close-packed planes,” because they minimize the 
spacing between atoms; they are the most common slip planes because the atoms do not 


Table 2-1 Atomic density on FCC crystal planes. 


FCC Plane 

Atom/Area 

Atom/Area 

{100} 

2/a 0 2 

2 

{110} 

2/ V2 a 0 2 

1.414 

{111} 

4/ i/3 a/ 

2.309 


Source: Dieter, 1986. 


have to travel great distances to reach another atomic position. Close-packed directions, 
like close-packed planes, minimize the distance between atoms. Therefore slip often 
occurs along close-packed directions in close-packed planes to minimize the amount of 
energy needed for displacement. 

Elasticity Matrix 

The energy needed to move in any direction is related to the elastic constants, and 
any material can be completely defined with 36 separate elastic constants. However, 
most material structures obey some type of symmetry, which reduces the number of 
independent constants. As indicated in Table 2-2 isotropic materials have only two 




Table 2-2 Symmetry in various crystal structures 


Crystal 

Structure 

Rotational 

Symmetry 

No. Independent 
Elastic Constants 

Tetragonal 

1 fourfold rotation 

6 

Hexagonal 

1 sixfold rotation 

5 

Cubic 

4 threefold rotations 

3 

Isotropic 


2 


Source: Dieter, 1986. 


independent elastic constants. Two of the three: E, modulus of elasticity; G, shear 
modulus or o, Poisson’s ratio will completely define the material properties in any 
direction. Cubic structures are highly symmetrical and require only one additional elastic 
constant, which reduces the overall elasticity matrix considerably. 

The process of reducing the original 36 constants down to three independent 
cubic constants begins with the original full elasticity matrix, [ay]: 

a H ai 2 ai 3 a 14 a 15 a 16 
a 21 a 22 a 23 a 24 a 25 a 26 
r , a 31 a 32 a 33 a 34 a 35 a 36 

[ay] = (2-4) 

a 41 a 42 a 43 a 44 a 45 a 46 

a 51 a 52 a 53 a 54 a 55 a 56 

a 61 a 62 a 63 a 64 a 65 a 66 

The subscripts correlate to the stress and strain components; the strain then is: 

ejsajjOj (2-5) 

CHAPTER 3Assuming elastic potential exists (i.e. isothermal deformation) the following 

relationship is achieved in equilibrium: 


Therefore the most general matrix is reduced to 21 components: 




a ll a 12 a 13 a 14 a 15 a 16 

a 12 a 22 a 23 a 24 a 25 a 26 

a 13 a 23 a 33 a 34 a 35 a 36 

[ay] = (2-7) 

a 14 a 24 a 34 a 44 a 45 a 46 

a 15 a 25 a 35 a 45 a 55 a 56 

a 16 a 26 a 36 a 46 a 55 a 66 


CHAPTER 4Nickel-base superalloys are orthotropic, meaning they have three orthogonal 

planes of elastic symmetry. Including the effects of cubic elastic symmetry, also called 

“cubic syngony” (Lekhnitskii, 1963) reduces the number of independent constants to 

three for the final elasticity matrix: 

a ll a 12 a 12 0 0 0 

a 12 a ll a 12 0 0 0 

a 12 a 12 a l 1 9 0 0 

0 0 0 a 44 0 0 

0 0 0 0 a 44 0 

0 0 0 0 0 a 44 


[ay] = 


( 2 - 8 ) 


The constants are defined by the modulus of elasticity, the shear modulus, and Poisson’s 
ratio along given directions: 



(2-9) 


1 

a 44 = — 
Gyz 


( 2 - 10 ) 


Vyx _ v xy 
Fxx Pyy 


( 2 - 11 ) 


For Material A the reported elasticity constant values are (Swanson and Arakere, 2000): 


a n = 6.494 E - 8 


344 = 6.369 E — 8 


aj 2 = -2.603 E - 8 


(psi) 


Using these elastic constants and the direction cosines for a given orientation, it is 
possible to calculate the modulus of elasticity of the material in any direction (Dieter, 
1986): 


1 W \ ] 2 2 2 2 22 

= ai | - 2- (a ! ! — a 12 ) - - a 44 • a 3 |3 3 + P 3 y 3 + a 3 Y3 

H1VW 2 


( 2 - 12 ) 


Where a 3 , p 3 and y 3 are the direction cosines from the load direction to the x, y and z- 
axes, respectively (often called /, m and n ). Some materials, such as tungsten, are 
essentially isotropic, even in single-crystal form, because their elasticity is nearly 
constant in any direction. Others, like nickel-base superalloys, vary considerably. 


Notch Tip Deformation 

The behavior of an ideal material, one that has an infinite body free from any 
irregularities or stress concentrations, is very different from the typical superalloy 
applications. A turbine blade, for example, has a complex geometry, is exposed to 
mulitaxial, centrifugal, and contact stresses; and must withstand extreme temperature 
gradients, among other complex conditions. In an attempt to study the state of stress for 
more complex specimens, notched tensile specimens are often used to represent either 
areas of stress concentration or theoretical fracture conditions. 


Rice (1987) 

Rice provided the foundation for much recent and current work in the area of crack/notch 
tip stress and strain analysis by examining the mechanics of both FCC and BCC (body- 
centered cubic) notched specimens loaded in tension. Rice introduced the use of a plane 



strain assumption, which he combined with the CRSS criteria to predict specific zones, or 
“sectors.” These angular regions radiating from the notch tip characterize different 
regions where specific deformation mechanisms are at work. He began with two 
commonly studied orientations (later further studied by Crone and Shield, Figure 2-11) 
and derived an analytical solution to predict the active slip systems around a notch, as 



Figure 2-5 Notch direction terminology. 

well as the specific angles defining the boundaries between sectors. Rice, and others who 
follow his work, typically refer to tensile test specimen orientations in terms of notch 
plane, notch growth direction, and notch tip direction (Figure 2-5). His terms will be 
maintained throughout this discussion to prevent confusion. 

Based on the plane strain assumption, the yield criteria can be presented in terms 
of a “yield area,” which outlines the boundaries where the state of stress will correspond 
to the yield stress of the material (Figure 2-6). The slip activity is limited to certain 
allowable slip systems, which combine to produce large in plane strains in accordance 
with the plane strain assumption (i.e. £ zz = y yz = y zx = 0). The plastic field Rice constructs 



Figure 2-6 Yield surface based on plane strain component stresses. 

Source: Rice, 1987. 

at the crack tip also correlates to a plane strain state of stress in an isotropic material. By 
essentially eliminating the effect of single crystal anisotropy, this solution likely would 
not correlate well to experimental data. Another effect of the simplified elasticity of 
Rice’s solution is the inability to solve for a detailed strain field based on the state of 
stress near the tip. Assuming the out of plane stress and strain is equal to zero near the 
notch tip, the effect of the notch on creating a triaxial state of stress is essentially 
nullified. 

Rice’s solution represents a continuous solution in terms of the radius and angular 
displacement from the tip, where the state of stress is constant within each sector. The 
solution is symmetric about the growth axis for each case based on crystal lattice 
structure, so only the positive half-plane will be discussed. Rice notes the only way all 
the notch surface boundary conditions can be met for all angles is if 1) the stresses in 
some sectors are below yield, or 2) there are discontinuities at certain angles. Rice 
assumes a perfectly plastic stress field, so all sectors must be at yield or past, and the first 
condition cannot be changed. (Note: this assumption is contrary to varying elastic and 



plastic sectors commonly observed during fracture testing.) Therefore the boundaries are 
defined as the radial lines where a discontinuity occurs at specific angles where the slip 
shifts from one system to another. (Note: Rice also solves for the second condition, 
maintaining continuity in order to later apply his solution to a separate crack growth 
model. For this case, some angular sectors may be elastic while others are plastic, 
corresponding to those locations where the solution intersects the aforementioned yield 
surface.) 

Rice’s solution makes no distinction between the two orientations’ sector 
boundaries, or between FCC or BCC crystal structure. Both crystal orientations and 
structures predict boundaries at 55°, 90°, and 125°. The slip systems do change between 
orientations and corresponds to a switch between the slip system shear and normal 
directions; thus orientation predicts no effect on the yield surface or sector boundaries. 
Rice notes the weakness of this attribute, based on contradictory experimental studies, 
which is tied to the actual rotation of the crystal lattice he has ignored. He also 
acknowledges the simplification of the plane strain assumption, which was intended only 
as a starting point. Finally, Rice neglects strain hardening, but encourages incorporating 
all simplifications into future models. 

Shield. CuH 996) 

Shield conducted several tests of notched single crystal copper specimens to 
correlate Rice’s analytical models with his experimental observations. Based on 
symmetry and the previous work done by Rice, Shield chose to test two identical 
specimens with a notch in the [Oil] plane and a [100] notch direction (Figure 2-7, 
matching with the later “Orientation II” in work by Crone and Shield, 2001); the sectors 



and boundaries are symmetrical about the [100] axis, therefore only the upper half-plane 


(+0) will be discussed. Following Rice, four sectors were predicted with boundaries at 0°, 
55°, 90°, 125°, and 180°. Recall, Rice chose these angles based on his plane strain 
assumption: the active slip planes must actually represent a combination of slip systems 
that result in plane strain deformation for his assumption to be valid. Furthermore, the 



Figure 2-7 Specimen orientation for Shield’s Cu test specimen. 


sector boundaries can only occur at angles that correspond to the active slip system 
directions or their normals within the adjacent sectors. For example, if the active slip 
directions in a given sector are at 55° and 125° to the [100] axis, then that sector’s 
boundary angles must be from the set: 55°, 125° (both slip), -35°, 145°, 35° or 215° (all 
normals). One conclusion that must follow from Rice’s assumptions is that 0° and 180° 
cannot be actual boundaries unless the slip or normal matches those angles, which it does 
not (Figure 2-8). Therefore, the first and fourth sectors really span from -55° to 55° and 
125° to -125°, respectively, which gives six sectors total in both half-planes versus eight. 



Slip at 55° 

Slip at 125° 

Slip at 180° ..... 

Figure 2-8 Slip systems predicted by Rice. 

Source: Crone and Shield, 2001. 

The experiments done by Shield, and Crone and Shield, use four-point bending to 
deform the single crystal specimens plastically near the notch tip (Figure 2-9). Moire 
interferometry analysis is then used to determine the strain fields and sectors. Shield’s 
prediction of slip sectors, like Rice, is then based on plastic strain fields, rather than 
yield-inducing shear stresses. Shield assumes the elastic strain is negligible and the total 
strain is therefore equivalent to the plastic plane strain (i.e. a perfectly plastic state). The 
interferometry analysis technique uses a grid-like diffraction grating applied to the 
surface of the specimen. A laser works in combination with several mirrors to resolve 
the strains into u-field and v-field component fringe patterns (Figure 2-10), which are 
used to calculate the actual strain values. Moire interferometry is a well-established 



Figure 2-10 Moire interferometry strain field. Shield determined the sector boundaries based on 
the obvious visible patters; line emphasis mine. 

Source: Shield, 2002. 


Shield experimentally tested two separate specimens and measured the boundary 
angles, based on strain fields (Table 2-4). Applying several load levels, beginning at the 
yield point, he observed whether the load had any effect on the sector boundary angles. 


Table 2- 4 Experimental sector boundar y angles. 


Sector 

Boundary 

Angle 
0 (deg) 

1-2 

43 

2-3 

62 

3-4 

100 


No significant sector changes were observed as the load increased, indicating that the 
sectors are independent of the load level. However, as the load increased the amount of 
plasticity near the notch tip increased as well, allowing for some stress relaxation as the 
material responds. The high plasticity near the notch tip also has an impact on the 
method of measuring the sector angles. To avoid the large deformations near the tip, 
which do not represent the desired strain fields from which the sectors are determined, 
the origin for angular measurements are offset by 50pm and the strains are presented 
along a 400pm arc. Note: with the notch’s mean width of 100pm these measurements are 
taken at 8*p, where p is the notch radius. This is a significant distance from the notch 
tip, regardless of elastic versus plastic conditions, and seems disproportionate to get an 
accurate account of the behavior due to the notch. 

Shield’s results from the low load levels show similarities to Rice’s model, 
however Shield’s experimental results do not correlate to Rice’s model at high plastic 




strains (Figure 2-8 and Table 2-4). As would be expected in a plastic model, these strains 
are not proportional (i.e. they do not have a linear response). In fact, they do not even 
maintain the same order in their relative level of activity: Sector 2 initially shows the 
most activity, while at higher strains Sector 3 becomes more dominant. This changing 
slip activity (with load level) seems inconsistent with the constant sector boundaries. 
Constant angles indicate the boundaries are independent of the active slip within the 
sector; perhaps the boundaries appear constant due to relaxation in the material between 
the load steps. Between the two samples the boundary angles were similar but not 
constant, highlighting the variance that can result from the material structure alone and 
any flaws that may be present, despite a constant specimen orientation and test 
conditions. The precise notch geometry also a factor is, which is extremely difficult to 
duplicate accurately. 

Plastic deformation at large strains is easily observed through slip lines; Shield 
used this visible slip evidence to compare the strain sectors determined by Moire 
interferometry. The resulting sector boundary angles matched well with the strain field 
images. Since the activated slip systems are a function of the RSS, and because they 
correlate well to interferometry results, the sectors determined by strain field analysis 
appear to be analogous to those determined by a stress-field analysis. Shield’s 
experimental results, and their dissimilarity to Rice’s analysis, again highlight the need 
for a more accurate predictive model. 

Crone and Shield, Cu and Cu-Be (2001) 

Extending the work of Shield (1996), Crone and Shield continued experimental 
studies of notch tip deformation in two different orientations of single crystal copper and 



copper-beryllium tensile specimens (Figure 2-1 1). (Orientation II will be presented first 
due to the widespread use of this orientation in other studies). Again using Moire 
interferometry strain fields and visible slip evidence to determine sectors, sector 
boundaries are defined as the obvious visible changes from one area to another. The 
visible slip patterns determine slip activity, but as the authors note, a lack of visible slip 
does not rule out any activity. Slip systems may be activated internally, rather than at the 
surface, or may show varying patterns on the surface as deformation continues. 


[ 010 ] 



[ 101 ] 


[ 101 ] 


[ 010 ] 


Orientation I 


Orientation II 


Figure 2-1 1 Orientations experimentally tested by Crone and Shield. 


Crone and Shield compared their experimental results with Rice’s analytical 
solution, as well as numerical FEA solutions by Mohan, et al. (1992) and Cuitino and 
Ortiz (1996) (Table 2-5). Both numerical models run on the plane strain assumption, 
although Cuitino and Ortiz later conclude the problem cannot truly be plane strain 


Table 2-5: Orientation II sector boundary angle comparisons. 


Sector boundary 

Experimental 

Analytical 

Numerical 

‘Boundary angles 
in degrees 

Crone and Shield 
(2001) 

Rice 

(1987) 

Mohan, et al. 
(1992) 

Cuitino and Ortiz 
(1996) 

1-2 

50-54 

54.7 

40 

45 



2-3 

65-68 

90 

70 

60 

3-4 

83-89 

125.3 

112 


4-5 

105-110 


130 

135 

5-6 

150 





Source: Modified from Crone and Shield, 2001 



Figure 2-12 Experimental slip sectors from Crone and Shield. 
Source: Modified from Crone and Shield, 2001 . 


due to large strain differences internally and at the surface. Even with the same plane 
strain assumption as Rice, the numerical and analytical models do not match; all three 
differ from the experimental results (Figure 2-12 and Figure 2-13). Like Shield’s earlier 
results, the experimental results are somewhat ambiguous due to the “annulus of 
validity," where Crone and Shield take their measurements (Figure 2-12). This annulus, 
following Shield and Kim (1994) corresponds to the radial area from 350-750pm from 
the notch tip. The notch width is between 100-200pm, making the notch radius between 




50- 100pm. Therefore the annulus, and the region where the sectors are measured, is 
anywhere from 3. 5-7.0 and 7.5-15.0 times the notch radius from the tip. These distances 
would place the sectors well out of the range of any elastic deformation and clearly can 
only be used where extreme plastic deformation exists. However, Crone and Shield 
preserve this annulus to avoid material too close to the notch tip, which is dominated by 
geometry, and a larger annulus which would begin to impose on these particular 
specimens’ boundaries. 

Crone and Shield describe the observed slip activity as it begins in a single sector; 
as deformation proceeds, more slip lines become visible in the same sector at further. 


1 



Slip at 180° 


Figure 2-1 3 Full field slip sectors and slip line traces from Crone and Shield. Note: these sectors 

differ from Rice (1987) and Shield (1996). 

Source: Crone and Shield, 2001. 


radial distances from the notch. This result is intuitive because a larger load will allow 
higher stresses even as distances from the notch tip increase. The number of slip lines 
then grows until the sector is “filled”. Although it is unclear what exactly is meant by 
“filled”, it is clear that, either the initially visible slip system, or another with the same 



visible trace, remains active throughout the experiment. Crone and Shield clearly 
observe horizontal slip traces directly ahead of the notch (Sector 1) for Orientation II, 
however they discount their observations and instead label the slip as “elastic” in order to 
compare their solution to other perfectly plastic sharp crack solutions (Figure 2-13). 
Citing Saeedvafa and Rice (1989) they further explain these traces as a function of 
hardening and not plastic deformation, but they do not make such accounts for the other 
sectors. The only other sector determined to be elastic is Sector 5, and this is due to its 
lack of any visible slip activity as well as low strain measurements. Consequently, their 
determination of an elastic zone directly ahead of the tip is not consistent with their 
visible experimental results. As predicted, the experimental sectors are symmetric about 
the [101] axis. Only Sectors 3/-3 and 4/-4 show variation between +0 and -0 (Figure 
2-13); the apparent asymmetry in these sectors is likely caused by specific notch 
geometry/irregularities, material defects, etc., or by a slightly inaccurate notch plane 
direction. 

Orientation I represents Orientation II rotated by 90° about the notch tip direction, 
such that the notch plane and the notch tip/crack growth directions are switched. 

Contrary to the equivalent sectors predicted by Rice, Crone and Shield’s observed sectors 
show a marked difference with orientation, varying in both specific boundary angles and 
in the number of sectors (Table 2-6). Along with Rice’s prediction, none of the 
numerical solutions correspond to Orientation II, which they were designed for, or for 
Orientation I. The numerical finite element models, however, agree more closely than 


Ta ble 2-6 Sector boundary angles for Crone and Shie ld. 


Sector boundary 

Orientation 

1 

Orientation 

II 

‘Boundary angles 
in degrees 

(101) Plane 

(010) Plane 

1-2 

35-40 

50-54 

2-3 

54-59 

65-68 

3-4 

111-116 

83-89 

4-5 

138 

105-110 

5-6 


150 


Source: Crone and Shield, 2001 


the others. Despite the disagreement of the experimental work to analytical and 
numerical solutions, Crone and Shield again support a plane strain assumption based on 
the numerical results from Cuitino and Ortiz (1996). The FEA results on a central plane 
of the model appeared to correspond more closely to the experimental results. Crone and 
Shield assert this agreement to the central plane FEA results suggests that plane strain is 
accurate for specific locations. However, obtaining results only from a central plane does 
not represent the true specimen being modeled. In fact, if there is a preference it should 
be for the model to be more accurate at the surface; where slip is thought to originate and 
actually observed, allowing a better comparison. 

Schulson and Xu, NFA1 (1997) 

A recent study by Schulson and Xu (1997) examined the state of stress at a notch 
tip for single crystal N^Al, the y’-component of Material A. Using three-point bending, 
two crystals were deformed at first elastically, then plastically until a crack formed, and 
finally through a small degree of crack propagation. Both samples were oriented with the 
growth direction [010], the notch plane direction [-10-1] and the notch tip direction 




[- 10 - 1 ] 



Figure 2-14 Schulson and Xu specimen orientation. 

[-101] (Figure 2-14). To solve an analytical model to compare their experimental results 
to, they assumed elastic isotropic material conditions and calculated the stress field 
around the notch based on the equations for a sharp notch (Anderson, 1 995). Although 
these equations may have been accurate for previous experiments they conducted with 
polycrystalline NisAl (Xu and Schulson, 1996), the isotropic assumption is not valid for 
their single crystal specimens. Schulson and Xu carry out two sets of stress calculations, 
one based on plane stress conditions and another with plane strain. The plane stress 
assumption can be used for large, thin plates with in-plane loading, where only a x , a y , 
and x xy are present (Bickford, 1998); the thickness must be small compared to the height 
and width of the specimen. However, as Schulson and Xu note in their paper, the notch 
causes a triaxial state of stress, where out of plane material stresses exist. Therefore, the 
plane stress assumption is not valid for a single crystal material. The plane strain 
assumption also relies on isotropic material conditions and can be applied when the out 
of plane strain is zero. A typical application where the plane strain assumption is applied 
is a hollow cylinder, in which the radial strain is assumed to be zero when the length of 
the cylinder is sufficient. However, for anisotropic single crystal materials the strain is 



greatly affected by the material orientation. The simplified stress tensor matrix (Eq. 2-8) 
is valid only when the specimen is loaded parallel the FCC lattice edges (Swanson and 
Arakere, 2000). Assuming the direction of their test specimens is accurate, Schulson and 


Xu may use this basic transformation. However, the true component stresses will have 


out-of-plane components; therefore plane stress is not valid. Using Eq. 2-5 the strain is 


equal to the stress tensor matrix multiplied by the stress: 
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Since out of plane strain is present (e 2 , y 2X , y yz ) the plane strain assumption is not valid. 
Cuitino and Ortiz (1996) also came to this conclusion after observing differing stress 
fields at the interior and surface of their FEM, as well as slip that would produce out of 
plane strains. 

Using the plane stress and plane strain assumptions, Schulson and Xu calculated 
two sets of resolved shear stresses for the 12 primary “easy slip” systems. They then 
normalized these resolved shear stresses with respect to Ki, the mode 1 stress intensity 



factor, and r, the radial distance from the tip. Plotting the normalized shear stresses 
against the angular displacement from the notch tip resulted in a series of shifting 
maximum stresses, corresponding to the different sectors (Figure 2-15 and Figure 2-16). 
Although they achieved similar slip system results for each assumption (Table 2-7), the 
exact systems differ in certain sectors (II, III, IV, V), and the sector angles are very 
different in some cases (I/II, IV/V). The sector I slip systems are equal under both 
assumptions, but the stress is different. The difference in the specific slip systems in 


^ MO-1] 



Figure 2-15 Slip sectors under [-10-1] load. Actual experimental results. 
Source: Modified from Schulson and Xu, 1997. 


sectors II-III is of interest, but the systems under plane strain are merely the symmetrical 



systems from those of plane stress, except lib. In sector IV however, although the 
maximum stresses are resolved onto the same planes, the directions do not match 
between the two conditions. Finally, in sector V the plane strain stresses jump back to 


Table 2-7 Slip sectors for plane stress and plane strain assumptions. 


Sector 

Plane Stress 

Plane Strain 


6 (deg) 

slip systems 

0 (deg) 

slip systems 

1 

0-43 

(11-1)[-110] 

(-111)[01-1] 

0-23 

(1 1-1)[-1 10] 
(-1 1 1 )[01-1] 

II 

43-60 

(1 1 1)[1-10] 
(1 1 1)[01-1] 

23-60 

(-1-1-1)[-1 10] 
(-1-1-1)[01-1] 

III 

60-103 

(1-1 1)[1 10] 
(1-1 1)[01 1] 

60-107 

(-11-1)[-1 -10] 
(-1 1-1)[0-1-1] 

IV 

103-180 

(1 1-1)[01 1] 
(-1 1 1)[1 10] 

107-133 

(1 1-1)[101] 
(-1 1 1 )[101] 

V 

— 

— 

133-180 

(-1-1-1)[0-1 1] 


Source: Modified from Schulson and Xu, 1996. 




Figure 2-16 Maximum RSS slip system plots. Note the difference between the plane stress and 

plane strain assumptions. 

Source: Schulson and Xu, 1997. 


the {111} systems, as in sector II, yet at these same angles the plane stress is still on the 
{-111} planes (Table 2-7). The data correlates well for the elastic region, the “early slip” 
(sectors II and III) between approximately 43-100°, but greatly diverges in stress level in 




sectors I, I/II, and I V/V. Although Schulson and Xu use alternate assumptions of plane 
stress and plane strain for their component stress calculations under isotropic conditions, 
the close agreement between the two predicts initial slip on the { 1 1 1 } and {-111} planes. 


Table 2-8 Slip sectors for plane stress and experimental results. 


Sector 

Plane Stress 

Sector 

Experimental 


0 (deg) 

slip systems 


0 (deg) 

slip systems 

1 

0-43 

(1 1-1)[-1 10] 
(-111)[01-1] 

1 

0-38 

(11-1) 

II 

43-60 

(1 1 1)[1-10] 
(1 1 1)[01-1] 

II 

38-58 

(111) 

III 

60-103 

(1-1 1)[1 10] 
(1-1 1)[01 1] 

III 

58-100 

(1-11) 

IV 

i 

103-180 

(1 1-1)[01 1] 
(-1 1 1)[1 10] 

IV 

100-notch 

(11-1) 


Source: Modified from Schulson and Xu, 1996. 


Experimental results after significant plastic deformation reveal results that 
deviate from those predicted by either plane stress or plane strain, but are remarkably 
closer to the plane stress assumption (Table 2-8). Note, under the plane stress 
assumption, out-of-plane strain is still present. Sectors II and III underwent early slip 
nearly equally in both the +0 and -0 half-planes (which are symmetrical), with Sectors I 
and IV being activated at later stages. Although the slip plane family can be determined 
visually, some specific slip systems appear the same at the surface; in these cases 
transmission electron microscopy (TEM) must be employed to detect the specific slip 
system. Using TEM at the Will boundary, Schulson and Xu observed [01-1] and [Oil] 
dislocations, which confirms the good correlation to the plane stress model. 

Noting the scarcity of work in this field Shield (1996) calls for further 
experiments, particularly in pure tensile testing (versus three/four-point bending), to 




observe the details of the applied stress on the strain field. Also of interest is the 
observation of slip activity ahead of the notch, versus at large angles from the notch tip; 
numerical calculations predict plasticity in both regions, but Shield observed a higher 
degree of plasticity largely ahead of the notch tip. An ideal test would incorporate 
parameters for specimen size, type of test, plasticity, hardening, lattice rotation, etc., and 
eventually a fatigue crack rather than a notch. However, before incorporating such 
complexity the basic model must be better understood. Based on the assumption that 
stress and strain fields will correlate well, and on the assumption that yielding is a stress- 
controlled process alone (CRSS), this thesis builds off previous work by introducing an 
elastic analysis of a notched single crystal superalloy material. The resulting resolved 
shear stresses predict active slip systems and can be used to determine sector boundaries. 
The model presented here also branches into the tensile testing setup, verses the major 
bending models created thus far. 



CHAPTER 3 

ANALYTICAL PROCEDURE 


This study’s goal is to find the state of stress in the material coordinate system of 
a notched single crystal superalloy test specimen, and then use those stresses to calculate 
the resolved shear stresses in the 12 primary slip systems. For an isotropic material a 
single elastic constant governs the transformation from stress to strain (Eq. 2-5). 

However, the elasticity of an anisotropic material is not constant, and analyzing the stress 
and strain fields in a single crystal material must first begin with determining the specific 
material properties in a desired orientation. The stress-strain relation for an anisotropic 
solid with cubic symmetry has three independent constants in the material coordinate 
system, as well as a “stress tensor matrix,” instead of a single elasticity constant, that 
varies with orientation. Since test specimen orientation can vary, we need requisite 
equations to convert between the specimen and material coordinate systems. 

Coordinate System Transformation 

The first step in defining the elasticity matrix is to determine the precise 
orientation of the actual specimen, either in terms of the material Miller indices (direction 
indices) or angular measurements. In a physical material test specimen, it is nearly 
impossible to cut a sample such that the x, y, and z test axes are perfectly aligned to the 
material axes: [100], [010] and [001] respectively. Therefore, a transformation is needed 
to convert the known specimen stresses to the material coordinate system. The applied 




stresses are a measure of force per unit area, and are a product of the slip plane geometry. 
In other words, the stresses vary by slip plane and direction according to a particular 
plane’s cross-sectional area and its given orientation within a unit cube of the material. 
Therefore, these stresses are not affected by anisotropy. However, only two material 
properties are independent without anisotropic effects; the shear coupling induced in the 
three-dimensional model, and resulting component stresses, will not be properly 
accounted for without the third independent constant to define the single crystal material. 
Lekhnitskii’s book on the Theory of Elasticity of an Anisotropic Elastic Body is a good 
reference on this subject. 

The transformation equations presented here follow the procedures outlined in 
Lekhnitskii (1963) and Stouffer and Dame (1996). The transformation from the 
specimen to the material coordinate system can be accomplished by two methods. In the 
first approach, the angles between the original and transformed coordinate systems may 
be directly measured to find the direction cosines. This method is preferable if the angles 
are easily found. The second method, based on rigid body rotations, may be used for 
more complex orientations, where the angles between the two coordinate systems are not 
as obvious; here the Miller indices of the transformed axes must be known. In this 
process, the axes are rotated through a series of steps to arrive at the final transformed 
destination. As far as preference, neither method is more correct than the other, but in 
experimental specimens it is often quite difficult to determine the exact Miller indices; 
and it can be more convenient to measure the angles. 

Coordinate transformations may be performed as long as the orientation of the 
sample is known, and the transformation matrices can then be used to determine the 



stresses and strains resolved on any given plane or slip system. The original coordinate 
system will be defined as the material coordinate system: x 0 = [100], y 0 = [010], and 
z 0 = [001] (Figure 3-1). The transformed coordinate system is defined as the specimen 
coordinate system, and is offset by some angular displacement from the original axes. 
The specimen axes are denoted by x”, y” and z”; all properties associated with the offset 
system are also denoted by the double-prime symbol. 

Three-Step Coordinate Axes Transformation 

As noted above, the initial, or material, coordinate system is denoted by x Q , y 0 , and 
z 0 -axes: 



Figure 3-1 Material (x 0 , y 0 , Zo) and specimen (x”, y”, z”) coordinate systems. 

The easiest way to complete the total transformation to the specimen coordinate system is 
to break it into several rigid body rotations; here we will use a three-step process. The 
first transformation, to the x, y and z axes, is performed by rotating by vj/ 0 about the z 0 - 
axis (positive is defined as x 0 toward y 0 ) (Figure 3-2). 


y 0 [010] 





Writing the transformation for step one in matrix form is: 

X cos(v 0 ) sin(\|f 0 ) 0 Xo 

y = -sin(\|/ 0 ) cos(\j/ 0 ) 0 • Yo (3-4) 

z 0 0 1 Zo 

The second transformation, to the x’, y’ and z’ axes, is accomplished by reflecting the 
load vector onto the x-z plane and rotating \j/i about the y-axis (positive defined as z 
toward x) (Figure 3-4). 



Figure 3-4 Second rotation about the y-axis. 


Following the same process as step one, the second transformation in matrix form is: 

x' cos ( mm) 0 -sin(vi) x 

y' s 0 10 ■ y (3-5) 

z' sin(\|/]) 0 cos(\|/i) z 

The final transformation to the specimen coordinate system, or the x”, y” and z” axes, 
occurs by rotating the y’ and z’-axes by \j /2 about the x’-axis (positive defined as y’ 


toward z’) (Figure 3-5). 



Figure 3-5 Third rotation about the x’-axis. 


The third transformation in matrix form is: 


1 0 0 x . 

0 cos (^ 2 ) s * n ( x l / 2) • y' (3-6) 

0 -sin(\j/2) cos(\|/2) z ' 

The total transformation can then be calculated by multiplying the three individual step 

matrices together (Note: the first transformation becomes the last one multiplied): 



where 


x " «1 Pi Y! 

y M = «2 P2 Y 2 • yo 

z ” a 3 p 3 y 3 Zo 


(3-7) 


«1 

Pl Yl 

1 0 

0 

cos (i|/ 1 ) 0 

-sin(y ]) 

cos(yo) 

sin(yo) 

0 

a 2 

P2 Y2 

= 0 cos(\|/2) 

sin(v2) 

0 1 

0 

-sin(yo) 

cos(y 0 ) 

0 

a 3 

P3 Y3 

0 —sin (^ 2 ) 

cos(\|/ 2 ) 

sin(\j/i) 0 

cos (y 1 ) 

0 

0 

1 


The resulting values represent the cosines of the angles between the material and 
specimen coordinate system axes (Table 3-1). 


Table 3-1 Direction Cosines 



X 0 

y 0 

Z 0 

X" 

CXl 

Pi 

Yi 

y" 

a 2 

P 2 

Y 2 

z" 

a 3 

P3 

Y3 


Solving for the direction cosines, or the final transformation matrix: 


<*i Pi Yl 

«2 P2 Y2 = i 
<*3 P3 Y3 


cos (\|/])cos(\|/ 0 ) 

sin(v|/ 2 ) sir»(\|/ ]) cos(\|/ 0 ) - cos (\|/ 2 ) sin(\j/ 0 ) 
cos(v2)sin(\|/i)cos(vo) + sin(\|/ 2 )sin(vj/o) 


cos(\|/]) -sin(vo) 

sin(\|f 2 ) sin(\|/ ]) sin(\|/ 0 ) + cos(\|/ 2 ) cos(y 0 ) 
cos(\j/2) sin(\|/i) sin(\|/ 0 ) - sin(\|/ 2 )-cos(\|/ 0 ) 


-sin(\|/i) 
sin (\|/ 2) cos (\|/ ] ) 
c:os(vj/2) cos(vj) 


( 3 - 9 ) 


When fewer than three steps are used, the selection of primary and secondary rotation 


axes is somewhat arbitrary and it is important to verify the results. Several checks based 


on perpendicularity can be calculated to ensure a proper orthogonal coordinate 
transformation has been performed (Appendix A). 


Example Transformation 

To illustrate a two-step transformation we will consider a load applied in the 
[213] direction (Figure 3-6). This transformation is reduced to two steps by eliminating 
the first 



x 0 [100] 




Figure 3-6 Load in the [213] direction. 


y [010] 



x [100] 


Figure 3-7 Two step coordinate transformation. 


step (i.e. \j/ 0 = 0), and step two begins by reflecting the load vector onto the x-z plane. 
The reflection shows a triangle whose sides are the x and z Miller indices: x = 2, z = 3. 
The first angular translation, \|/j, is: 


2 

\\f i = at an — 


( 3 - 10 ) 


y/i = 33.69° 

In the same way, the second angle, \\ h, forms a triangle with the hypotenuse, h , of the first 

reflection and the y-translation: h = V~i 3 ? y = 1 . Therefore, the second angular translation 
is: 


\j/2 = atan 


1 


(3-11) 


yji = 15.50° 



(Note: vj /2 must be negative for the [213] orientation.) Solving for the direction cosines: 

«i Yi 10 0 cos(\|/j) 0 —sin (\|/ j) 

«2 P2 Y 2 = 0 cos ('t , 2 ) sin(y 2 ) 0 10 (3-12) 

a 3 p 3 y 3 0 -sin(v 2 ) cos(y 2 ) sin(vi) 0 cos(yi) 

a 1 Pi V 1 0.832 0 -0.555 

a 2 0 2 Y 2 = -0.148 0.964 - 0.222 

ot 3 p 3 y 3 0.535 0.267 0.802 

Following Appendix A all accuracy checks confirm a proper orthogonal transformation. 
Stress and Strain Transformation 

Once the direction cosines between the material and specimen coordinate axes are 
known, the load conditions can be applied and incorporated into separate matrices to 
correctly transform the individual stresses and strains. These transformed matrices are 
then used to solve for the resolved stresses and strains in each desired slip system. 
Following Lekhnitskii (1963) the stress transformation is: 

{a”} = [Q']{a} (3-13) 

{a}=[QT'{a”}=[Q]{o”} (3-14) 


Here [Q] is the stress transformation matrix, a function of the direction cosines: 



2 

«1 

2 

«2 

2 

a 3 

2 a 3 a 2 

2-aj 013 

2 a2 a l 


Pi 2 

P2 2 

P3 2 

2 P 3 P 2 

2 p rp 3 

2 P 2 Pi 


2 

Yl 

2 

Y2 

2 

Y3 

2'Y3'Y2 

2 Yl Y3 

2 Y2 Yl 

( 3 - 15 ) 

PlYl 

P 2 Y 2 

P3Y3 

(p 2 Y 3 + P 3 Y 2 ) 

(P, Y 3 + P 3 Yl) 

(PlY2 + P 2 Y 1 ) 


Yr a i 

Y2«2 

Y 3 CX3 

(Y2 a 3 + Y3 a 2) 

(Yi-a 3 + Y 3 «i) 

(Yi-a 2 + Y2 a i) 



oq Pi 012 P2 a 3 p 3 (a 2 p 3 + a 3 p 2 ) (arP 3 + a 3 Pi) (a r P 2 + a 2 Pi) 



{e”} = [Q’ c ]{e} ( 3 - 17 ) 

{e} = [Q'c] '{£”} = [Q t ]{e”} ( 3 - 18 ) 

The strain transformation matrix, [QJ, is also a function of the direction cosines, but 
differs from the stress matrix, [Q]: 


2 

ai 

2 

a 2 

2 

a 3 

a 3 a 2 

a r a 3 

a 2 a\ 


Pi 2 

P2 2 

P 3 2 

P3P2 

Pi P3 

P2P1 


[Q e]= y ,2 

2 

Y2 

2 

Y3 

Y3Y2 

Y1Y3 

Y2 Yl 

( 3 - 10 ) 


2pi Vi 2 P 2 Y 2 2p 3 y 3 (P 2 Y 3 + P3 Y 2 ) (Pl Y3 + P3 Yl) (P 1 Y 2 +P 2 Y 1 ) 

2Y 1 • ot 1 2Y2 a 2 2y 3 « 3 (y2 «3 + Y3 «2) (Yl «3 + Y3 «l) (Yl-a 2 + Y2 a i) 

2a 1 P 1 2a 2 p 2 2a 3 p 3 (a 2 P 3 + a 3 p 2 ) (aj P 3 + o^ P]) (cqP 2 + a 2 Pi) 



According to Hooke’s law, an isotropic material’s stress and strain for a uniaxial state of 
stress are related by: 

a = Ee (3-20) 

The generalized Hooke’s law for a homogeneous anisotropic body is: 

{a} = [A u ]{e} (3-21) 

where the [a,j] and [Ajj] matrices are the elastic constant matrices and: 

[Ajj] = [a.j]-' (3-22) 

and [ay] is a symmetric matrix such that: 

[a,j] = [a„] (3-23) 

Therefore 


{e} = [a,j] {a} (3-24) 


and 


{£”} = [a’.j] {a”} (3-25) 


The elasticity matrix also undergoes transformation, but it remains symmetric. A 
maximum of 21 individual constants may be present, depending on the specimen 
orientation: 

[a’lj] = [Q] T [a,j][Q] (3-26) 

Now, once the component stresses are known in the specimen coordinate system, the 
above equations can be applied to solve for the component stresses in the material 


coordinate system. 


Slip System Shear Stresses and Strains 

The component stresses define the complete state of stress of the material, but 
these values alone reveal little about individual slip systems or RSS activity. The 12 
primary slip systems (Table 1-1) are defined by both a slip plane and direction. 
Incorporating their geometry is analogous to performing a separate transformation, and 
allows all individual resolved shear stresses to be determined. 


Resolved Shear Components 

Following Stouffer and Dame (1996), the resolution of stress to the slip systems is 
calculated: 

{x}=c[S]{o} (3-27) 

Here c and [S] are constants defined by the slip plane and direction: 


J h'j 2 + k’i 2 + l'rju'i 2 + v'j 


2 ,2 

+ VVj 


(3-28) 


Sj = ( h’j- u’i k'j-v'j l'j. w'f — w'j -v’j -u’j) 


(3-29) 


Here, [u’ v’ w’] is the slip direction and (h ’k’ P) is the slip plane (recall Figure 2-2). For 
the 12 primary systems, c is constant and the full RSS matrix becomes: 
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Solving Eq. 3-30 for the 12 primary systems: 
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(3-31) 


The shear strains are calculated in the same way as the stresses outlined above: 


{y} = c[S]{e} 


(3-32) 


The stress and strain fields, fully resolved onto the primary octahedral slip systems, are 
now known and can be used to predict slip within a particular system. 


CHAPTER 4 

NUMERICAL SOLUTION: FINITE ELEMENT METHOD 
Finite Element Model 

Finite element software can be used to model the specific geometries and 
orientations of the tensile test specimens of interest (Figure 4-1). Using the ANSYS 
finite element software (Version 5.7), three different samples were modeled in order to 
predict slip activity and sectors around the notch. Two of the three orientations are based 
on previous work (Rice, 1987; Schulson and Xu, 1997; Crone and Shield, 2001) and all 
three correlate to collaborative work between the Mechanical Engineering and Materials 
Science Engineering (MSE) departments of the University of Florida (UF). As a stress- 
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Figure 4-1 Finite element analysis specimens and orientations. 


based process, slip deformation should be predicted by the numerical model’s highest 
individual resolved shear stresses. The slip systems that are represented by the highest 


resolved shear stresses should then be observable as slip lines in the experimental test 
samples. 

A solid tensile specimen was first modeled with the given load conditions and 
then compared to the analytical solution outlined in Chapter 3, using MathCad 2000 
Professional, to verify the numerical procedure and results. Two notches were then 
incorporated into the model to represent the experimental test specimens. The finite FEA 
component stresses were taken from the material coordinate system, around the notch, 
and then used in the transformation equations to calculate the individual resolved shear 
stresses. Data was analyzed over a wide range of radial and angular distances to create a 
complete stress field, later used to draw conclusions on sectors and slip activation. 

Verifying the Finite Element Model 

The complete tensile specimen, without any notches, can be analyzed according 
to the procedure outlined in Chapter 3 to verify the initial finite element model. 

Although any model size should result in a correct analytical solution, the same 
dimensions as those used for the actual specimens (minus notch geometry) were used for 
consistency. The applied load used for all numerical models is 100 lbs (7028 psi with the 
specimen cross-sectional area). This load will result in stresses well below the yield 
point of the material to ensure a purely elastic model. Since the stresses vary linearly 
with the load, the actual load is not important as long as the material remains in the 
elastic regime. The initial check gives excellent agreement between for the component 
stresses (Table 4-1); 


Table 4-1 Analytical and numerical component stresses for Specimen A. 


O (psi) 

o x 

o Y 

o z 

°xv 

<7 VZ 

Oxz 

Analytical 

0 

0 

7027.9 

0 

0 

0 




Numerical 

3.6642 

3.6643 

7028.4 

-3.7176 

4.86E-03 

-2.49E-02 


% Error 



0.007% 




the percent error is well within acceptabl 

e limits, so the initial model a: 

Ffirms an accurate 


coordinate and stress transformation. Although the stresses attest to an accurate model, a 
better test of a correct anisotropic model is the component strains, which take the 
transformed stress tensor matrix into account (Table 4-2). The component strains in the 


Table 4-2 Analytical and numerical component strains for Specimen A. 


e 

Ex 

Ev 

e 2 

Exv 

Evz 

Exz 

Analytical 

-1.83E-04 

-1.83E-04 

4.56E-04 

0.00E+00 

0.00E+00 

0.00E+00 

Numerical 

-1.83E-04 

-1.83E-04 

4.56E-04 

-2.37E-07 

3.09E-10 

-1.59E-09 

% Error 

0.016% 

0.016% 

0.044% 





material coordinate system again show nearly zero error, affirming a correct numerical 
model. Each orientation was validated, using the same procedure outlined above, to 
confirm correct stress and strain transformations before proceeding with the notched 
model. Incorporating the notch geometry into the existing model will now complete the 
actual specimen model. 

Specimen Orientation and Geometry 

The three numerical models all utilize the same geometry, including simplified 
notch geometry, to observe the effects of orientation without other defect/size 
considerations. The finite element model (FEM) is limited to the body of the specimen, 
rather than the entire geometry including the end grips. The reason for the abridged 
model is two-fold: First, the mechanics at the grips are considerably different from the 
center of the tensile specimen and include the effects of tensile rig contact pressure, 
loading rate, etc.. In the experimental model, different deformation mechanisms are seen 
at the grips, and the sample can even fracture there first. The goal here, however, is to 




model those specimens which fail at the central area of the specimen and analyze those 
stresses of interest. Second, in an attempt to gain more accurate results grips may be 
changed or updated, as is the case currently in the MSE department. The numerical 
model will not be subject to the type of grip used or any variations it may cause. 

The geometry for all three models is based on the experimental counterpart to 
Specimen A, and was modified to simplify the geometry. Table 4-3 shows the actual 
specimen geometry of Specimen A, as well as the modifications made for the simplified 
FEM (Figure 4-2). The notch is modeled as the combination of a perfect rectangle and a 
perfect semi-circle. In reality the notch will likely have some angular offset with the 
horizontal, as well as some y-displacement offset from the specimen’s center. Also, the 
actual notch tip is not semicircular, but rather a smaller arc. In the numerical model, 



both notch lengths and heights were set equal to those of the largest actual dimension. 
Additionally, the notch radius was set equal to the notch height (i.e. one-half the width of 



the notch) to form a perfect semi-circle. These geometrical simplifications can be easily 
removed for more specific test results, rather than our general focus on orientation. 


Table 4-3 Actual (Specimen A) and finite element specimen geometry. 


Specimen Geometry 

(mm) 

Actual 

FEM 

Width 

5.100 

5.100 

Height 

19.000 

19.000 

Thickness 

1.800 

1.800 

Right Notch Length 

1.300 

1.550 

Left Notch Length 

1.550 

1.550 

Right Notch Height 

0.113 

0.113 

Left Notch Height 

0.111 

0.113 

Right Notch Radius 

0.045 

0.113 

Left Notch Radius 

0.055 

0.113 


Model Characteristics 

Material Properties 

The finite element model is a linear-elastic, orthotropic model. The ANSYS finite 
element software has several three-dimensional elements available to account for 
anisotropic or orthotropic material properties. Choosing these elements, and defining 
either the three necessary independent stress tensors (an, an, Zaa) or the three 
independent directional properties (E, G, and v), will accurately model a single crystal 
material. To calculate the stresses in any direction, the material properties must first be 
defined within the material coordinate system. In the FEM, the model is created around 
the global specimen coordinate system (Figure 4-3). Therefore, the proper direction 
cosines must be used to create the material coordinate system. ANSYS aligns material 
properties with the element coordinate system; therefore the element coordinate system is 



aligned with the material coordinate system so that the directional material properties are 


suitably applied. 


Applied load 



Figure 4-3 Global and material coordinate systems. The specimen is created around the 
global system (x 0 , y 0 , z 0 ) and the material system (x”, y”, z”) is later specified. 

Elements and Meshing 

The ANSYS elements chosen for the FEM are PLANE2 and SOLID95. After the 
three-dimensional solid model is created, the front face is meshed with the PLANE2 
elements (Figure 4-4). This front face has precise element sizing along the defined radial 
lines around the notch tip at 5° intervals. The PLANE2 element is a two-dimensional, 
six-node triangular structural solid. It models irregular structures well and uses accurate 
quadratic displacement functions (ANSYS 5.7 Elements Reference, 1999). The nodes 
also incorporate orthotropic material properties, however these nodes are later deleted 
after serving their sizing function and their material properties are not relevant. 


Figure 4-4 Finite element mesh about the notch tip. 


Once the front face is meshed, three-dimensional elements are swept through the 
volume to complete the meshing of the model. Working in conjunction with the two- 
dimensional elements on the front face, the three-dimensional elements retain their sizing 
definitions on each of the x-y planes through the specimen’s thickness (Figure 4-5). 
SOLID95 is a three-dimensional structural solid with 20 nodes; each node has three 
degrees of freedom with translation in the x, y, and z directions (Figure 4-6). It has the 
same basic structure as the ANSYS SOLID45, but the addition of mid-side nodes enables 
more accurate solutions near irregular shapes. SOLID95 is preferentially selected for this 
very reason, since the area around the notch is of the most interest. Also critical is the 
inclusion of orthotropic material properties, which again correspond to the element 
coordinate system. The original element structure is left intact, rather than enabling 



pyramidal shapes to reduce the solution run time. Although the pyramidal option is less 
costly it may result in less accurate solutions, especially where there are large stress 



Figure 4-5 Element sizing on the front face and through the thickness. The two dimensional 
frontal elements are swept through the thickness of the sample to retain guidelines for the three- 
dimensional elements. Left, showing the edge of the specimen; right, showing the elements 

through the thickness. 


gradients. SOLID95 also supports plasticity, creep, large deflections and large strains 
among others, and is well suited for further future development of the numerical model. 




Figure 4-6 ANSYS SOLID95 element. 


Source: ANSYS 5.7 Elements Reference, 1999. 


As described previously, ANSYS aligns the material properties with the elements, 
which the model aligns with the material coordinate system. ANSYS then transforms the 
stress tensor matrix to any desired direction for stress or strain calculations. Due to a 
weighted nature in the nodal stress solutions, stress calculations are only available at the 
corner nodes, rather than all corner and mid-side nodes. To observe both the near-field 
state of stress and those stresses at larger distances from the notch (radial and angular), 
six concentric arcs were created at the following radii from the notch tip: 0.25*p, 0.50*p, 
1.00*p, 2.00*p, 3.00*p, and 5.00*p; where p is the notch radius (Figure 4-7 and Figure 
1-5). (Note: Initially the 3.00*p arc was not included and a 10.0*p arc was present. 



Figure 4-7 Radial arcs used for element location and sizing; centered at the notch tip. 

However, the larger arc intersected with its counterpart around the opposing notch. The 
nodes selected at 3.00*p were not created along a defined radius, but were carefully 
chosen by hand through the ANSYS Graphical User Interface (GUI).) 





Solution Location 


The element sizing of the FEM allows data to be collected on any of five separate 
x-y planes (Figure 4-5), including the front, middle, and back planes. In the interests of 
our collaboration with the MSE department, the current focus is on the front surface of 
the specimen. Here, at the surface, we can observe slip lines and make comparisons to 
the numerical model. Unlike models that use a plane strain assumption, it is not critical 
to collect data on the central planes versus another location. The decision to examine the 
state of stress at the surface is practical both from a mechanics perspective (considering 
surface effects) and the desire for collaboration. 

The nodes that lie on the concentric arcs around the notch tip are selected by first, 
manually viewing the arcs in the GUI and then, choosing the corresponding nodes from 
0° to 180° from the x-axis (toward the y-axis). Calculations do not need to be performed 
from 0° to -180°, because the superalloy, an FCC material, exhibits symmetry about the 
<100> and <1 10> axes (see Chapter 2). To check this assumption, stress calculations 
were performed on Specimen A for the 0.25*p arc and confirmed the symmetry. For the 
negative angles, the particular resolved shear stresses change direction, along with the 
changing direction of shear deformation between the top and lower halves of the test 
specimen. However, the maximum values of the resolved shear stresses are consistent 
between the positive and negative angles, as are the planes with those stresses. 

Assumptions 

A linear elastic model may predict the initial slip, but cannot precisely predict 
subsequent behavior due to plasticity around the notch tip. Some other modeling 
assumptions include low temperature deformation, no account of microstructural (i.e. 



dislocation) mechanisms, and no crystal lattice rotation. Finite element analysis is 
capable of accounting for changing temperature effects. However, in order to simplify 
the model and collaborate more closely with the MSE department, the FEA applied 
material properties at a constant (room) temperature. The exclusion of microstructure 
has also been discussed. The last assumption, however, is valid for this elastic model; 
crystal lattice rotation generally requires a few percent strain in a single slip system to 
occur (Stouffer and Dame, 1996) and the load used here is extremely low. The FEA here 
assumes only elastic deformation, and seeks only elastic stress and strain trends rather 
than specific values at fracture. Future research can easily build from the current model 
to incorporate plasticity, as well as creep or other desired behaviors. 



CHAPTER 5 

RESULTS AND DISCUSSION 


The main purpose of using FEM was to determine the slip activation and resulting 
sectors for a notched single crystal tensile specimen. Results are presented here for an 
elastic stress response to the given load on the primary octahedral slip systems about a 
notch tip, with the RSS as a function of radial and angular position. The absolute values 
of the resolved shear stresses are used for easy comparison, since the direction is 
irrelevant for this analysis. As predicted, each of the three orientations shows different 
slip activation, and at different locations. The determination of sectors is a more complex 
question than initially indicated by the majority of the literature, and will first be defined 
as the dominant slip system at any theta (i.e. the system with the highest RSS at any 
theta). Each orientation will be presented individually at first and then compared with 
each other. 


Specimen A 

Recall, Specimen A was tested with a [001] load, a [-1 10] growth direction, and a 
[1 10] notch plane (Figure 4-1). The maximum RSS at any location is 12=25,000 psi at 
r = 0.5*p and 105°, or 3.56*a ap piied. Results are presented for the 12 primary RSS values 
from 0.25*p to 5*p and from 0° to the top of the notch (100° for 0.25*p up to 170° for 
5*p) (Figures 5-1 to 5-6). The slip system with the maximum RSS varies with radial and 
angular position; sectors were determined for each radius by the overall maximum RSS, 


68 



or the dominant slip system (Table 5-1). Notice that for a given angle the dominant 
system (and often the other activated systems) is not constant over the range of radii. As 
the state of stress changes away from the notch tip, the RSS also changes; however, each 
component stress does not change uniformly with the others, so their status relative to 
each other may change and shift with the magnitude of their stresses. Therefore the 
dominant slip system is likely to change one or more times for a given theta. 

When the RSS is plotted against theta the effect of theta is clearly shown. The 
effect of radius, on the other hand, can be more clearly seen by combining the results in a 
single plot with the 12 primary stresses for the entire range of radii (Figure 5-7). The 
stresses at each radius have been scaled so they appear in ascending order (with respect to 
radial distance) from the origin. The amount of data can be difficult to interpret, but 
shows the overall changes well; for more detail the individual radial figures should be 
investigated. To observe the effect of radius at a given angle, let us look at the RSS 
changes along 0 = 85° as an example: at r - 0.25*p, x 6 is the maximum RSS; however the 
largest stress quickly shifts to X2 from r = 0.5*p to 3*p, and then back to X6 at r = 5*p. 
Some angles do maintain a single dominant slip system for all radii (i2 at 57°-68° here), 
but overall the slip systems and sectors are highly variable throughout the RSS field. 
Overall the RSS field for Specimen A is dominated by ij, X2, and X3 on the (1 1 1) plane, 
and to a lesser extent X6 and xn, both on {-1 1-1 } planes. 
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Specimen B 


Specimen B was tested with a [-1 10] load, a [1 10] growth direction, and a [001] 
notch plane (Figure 4-1). The maximum RSS at any location is i4=25,720psi at r=0.5*p 
and 105°, or 3.66*o ap piied- Note, this stress is moderately higher than the maximum for 
Specimen A and is in a different slip system, however it does occur at the same location. 
Results are again presented for the 12 primary RSS values from 0.25*p to 5*p and from 
0° to the top of the notch. Like Specimen A, the slip system with the maximum RSS 
varies with radial and angular position. The dominant slip system sectors at each radius 
show a larger number of sectors than Specimen A, especially for higher angles (Table 5- 
2). Again, the dominant system is not constant over the range of radii. 

Like Specimen A, the resolved shear stresses change values and shift positions 
relative to each other with respect to theta. Combining and scaling the 12 primary 
stresses for the entire range of radii gives the entire RSS field (Figure 5-14). Although 
the number of activated slip systems is higher for Specimen B, the angular range that 
maintains a single dominant slip system for all radii is significant with 14 ranging from 
63°-100\ Overall the RSS field is dominated by 14, x 6 , t 7, and i 8 on the {-1 1-1} planes. 
The field also shows 19 and T12 briefly, again on the {-11-1} planes, as well as xi and T3, 
which emerge at high angles on the (1 1 1) plane. 
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Specimen C 


Specimen C was tested with a [1 10] load, a [001] growth direction, and a [-1 10] 
notch plane (Figure 4-1). The maximum RSS at any location is t2=22,200psi at r = 0.5*p 
and 105°, or 3.16*a ap piied- This stress is lower than the maximum for either Specimen A 
or B, again occurring in a different slip system but at the same location in the RSS field. 
Results are again presented for the 12 primary RSS values from 0.25*p to 5*p and from 
0° to the top of the notch. Like the two previous specimens, the slip system with the 
maximum RSS varies with radial and angular position. The dominant slip system sectors 
for each radius reveal a larger number of sectors than Specimen A at higher angles, but 
approximately the same as Specimen B (Table 5-3). Again, the dominant system is not 
constant over the range of radii. 

Consistent with Specimens A and B, the resolved shear stresses change values 
and shift position relative to each other with respect to theta. Combining and scaling the 
12 primary stresses for the entire range of radii gives the entire RSS field (Figure 5-21). 
Like Specimen B, although the number of activated slip systems is higher, the angular 
range that maintains a single dominant slip system for all radii is considerable with 12 
ranging from 64°-97°. Overall, the RSS field is dominated by xi, X 4 , t 5 , xio, and tjj on 
both the (1 1 1) and {-11-1} planes. At high angles X 3 , x 7 , and xg briefly appear, again on 
the (111) and {-11-1} planes. Specimen C, therefore, shows a dominant system on each 
of the four possible primary slip planes in the RSS field. 
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Table 5-3 Specimen C dominant slip system sectors. 
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Specimen/Orientation Comparison 


For each of the three specimens, the [001] orientation occupies a different axis of 
the notch geometry: for Specimen A [001] is the load axis, for Specimen B it is the notch 
plane axis, and for Specimen C it is the growth axis. Comparing these three orientations 
is advantageous, because although [001] is the most common growth direction (and 
presumably load direction) a complex part like a turbine blade may have stress 
concentrations close to one of the specimen notch orientations. 

The radius with the maximum RSS, r = 0.5*p, is a good location to begin to 
compare data. For all three orientations the number of dominant slip systems is very 
small at this radius: two for Specimen A and one for Specimens B and C. As noted 
above, Specimen B has the highest RSS overall, and it is also the only specimen with the 
maximum RSS not on the (1 1 1 ) plane. Of greater interest though, in terms of application, 
is that Specimen C has the lowest RSS. Not only does Specimen C have the overall 
lowest maximum RSS, but its second highest stress is significantly lower than the other 
two orientations as well. The “second maximum” for r = 0.5*p again occurs at 105° for 
all three orientations: For Specimen A 16=23,800 psi (3.39*a ap piied) J for Specimen B 
xi2=24,570 psi (3.50*a ap p| ied ), and for Specimen C x 5 =l 9,930 psi (2.84*a a p p iied). 

Recall, Specimen C is the only specimen to eventually activate all four primary 
slip planes, though at its maximum it is still limited to only the (1 1 1) plane. From a 
mechanics perspective, it is unclear whether the slip plane itself has an effect on the 
desirability (or lack thereof) of one orientation over another. In fact, certain attributes 
(such as the second, third, and fourth maximum RSS values) may actually be a more 
desirable criteria as they determine whether mechanisms such as cross-slip will occur. 
These mechanisms, though they may result in more deformation, can also avoid fracture 


by releasing energy through ductile deformation. Therefore, any propositions regarding 
the best orientation will be avoided here and left to those with a greater knowledge of 
dislocation and other atomic mechanisms to determine the most favorable design 
orientation. Nonetheless, from a merely stress-based approach Specimen C is the clear 
candidate to design for a larger load tolerance of the notched specimen. 

Experimental Results 

As noted in the introduction, these experiments had a dual purpose of analyzing 
the orientations tested, and also providing a numerical analysis to compare experimental 
results to, specifically those obtained by the MSE department. As of the writing of this 
paper, the MSE department had tested Specimen A in a paper published by Forero and 
Ebrahimi (TMS Annual Meeting, 2002); those results will be compared here. Forero and 
Ebrahimi applied a 1 100 lb tensile load to the specimen (exact geometry given in 
Chapter 2) whose orientation is the same as Specimen A. Although the specimen 
orientation is symmetric about the notch, some results vary between positive and negative 
theta (probably due to irregularities in the notch cutout); therefore all values are presented 
here. The 5*p radius will be used to minimize the effects of plasticity at the notch tip and 
correlate better with the elastic model. After scaling the numerical curves for the actual 
applied load, we can draw in the yield stress (x y i e id ~ 47 ksi) to see which slip systems are 
predicted to be activated (Figure 5-22). The FEA results predict slip activation from 0° to 
135°, however, in the experimental test activation is only visible at the given radius up to 
1 10°. Also, whereas several systems are predicted to be active in many locations, the 
experimental specimen generally shows a single dominant system in each sector (Figure 
5-23). Recall, the exact slip system cannot be determined by the slip trace analysis alone 



(Figure 5-24); nonetheless, correlating the known slip plane to the numerical prediction is 
still a good measure of the model’s accuracy. The single dominant slip systems predicted 
for each sector by the FEA and those indicated by the experimental specimen do not 
match (Table 5-4); however, the visible systems are among those predicted to be 
activated by the applied stress level. 



Figure 5-22 Experimental load for Specimen A. The dashed line indicates the yield stress of the 
material; any RSS curves above this line represent slip systems that should be activated. 



Table 5-4 Specimen A experimental results. 
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Figure 5-24 Tensile test specimen and surface slip lines. 
Source: Courtesy, Dr. Fereshteh Ebrahimi, UF. 


Application 

As shown in the previous section, a set of charts giving the RSS as a function of 
both 0 and p about a notch tip can be used to predict the active slip systems and sectors 
for a given load and yield stress. Because the resolved shear stresses are linear, the charts 
can be modified for any desired load level by simply multiplying the original applied load 
(100 lb = 7028 psi) by some constant to reach the new load level. The sectors can be 
determined by drawing in the yield stress level and obtaining those angles where the 
resolved shear stresses intersect that yield level. Any resolved shear stresses above the 
“yield line” should be activated and the intersections determine the sectors. 



Example 

A tensile test is conducted on a specimen with geometry and orientation identical to 
Specimen A. A 2 kip load is applied to the specimen, whose material is known to have a 
yield stress near 50 ksi. The original load, 1 001b ? is multiplied by 20 to reach the new 
load; an identical cross-section allows a similar multiplication of the originally applied 
7028 psi to 140.6 ksi. The original chart can then be scaled by the same factor, giving a 
range from 0 to 450 ksi. The yield line is drawn in at 50 ksi (Figure 5-25); all individual 
RSS lines above the yield should be activated at the corresponding intersecting angles. 
This example shows eight separate slip systems being activated for the entire range of 
angles at the given radius (Figure 5-25). Two additional slip systems are initiated at 40° 
and 42°, giving essentially two sectors with different possible activations: 0°-40° and 40°- 



Figure 5-25 Numerical plot adjusted for example load; yield stress indicated. 
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Figure 5-26 Numerical plot adjusted for example load; maximum RSS changes with theta. 


100°. However, the maximum RSS varies along the radial line, and the dominant slip 
systems shifts from in (-1 -1 1 ) to T 2 (1 1 1 ) to x 6 (-1 1-1 ), giving three sectors where the 
dominant system shifts: 0°-17°, 17°-82°, and 82M00°. Although the total number of 
activated systems and their sector angles is subject to the specific applied load, the 
systems with the highest RSS should always be the ones that undergo visible deformation 
first. Beginning at high angles, i6 should initiate first, followed by x 2 surrounding most 
of the notch, and finally ii i should initiate directly ahead of the notch. 

General Material Slip Activation 

Material A is a good experimental material to compare a simplified model to 
because it does not experience a large amount of strain hardening. Strain hardening, the 
result of cross-slip, and other deformation behaviors can inhibit analysis because the 



evidence of the original dislocation motion may no longer be evident. This case is well 
illustrated with another superalloy, “Material B,” which undergoes a large amount of 
cross-slip; cross-slip dominates the surface and makes it difficult to determine the 
original slip activation. The model’s elastic simplification also better predicts slip 
activation for low strain hardened materials because these materials do not undergo a 
large stage of plastic deformation before they fail. Therefore, the current model should 
correlate well with other materials with similar behavior. 

If the model is developed into an elastic-plastic model, then it should predict slip 
activation even more accurately in materials like Material A, and may even be applied to 
strain-hardened materials. However, since strain hardening is a result of dislocation 
mechanisms, whose scope exceeds the standard mechanics regime, it is unknown whether 
a model that does incorporate plasticity would accurately predict slip activation and 


sectors. 


CHAPTER 6 
CONCLUSIONS 

1. A three-dimensional linear elastic finite element model that includes the effect of 
material anisotropy provides a good indication of active slip planes and sectors at 
the surface of a notched tensile test specimen. 

2. Slip sectors determined by the stress field are not constant for a given material, as 
generally indicated by the literature, but are determined by the specific applied 
load and yield stress of the material. 

3. Under tensile loading, a specimen with a [1 10] notch plane (primary) orientation 
and a [001] notch growth direction will have lower resolved shear stresses than 
one with [00 !]/[-! 10] or [-1 10]/[1 10] notch plane/notch growth directions, 


respectively. 



CHAPTER 7 

RECOMMENDATIONS FOR FUTURE WORK 


The elastic model presented here appears to predict slip accurately, and can be 
used to predict fatigue behavior based on equations that incorporate individual RSS 
values, such as fatigue life (Swanson and Arakere, 2000). The RSS values, even though 
based on an elastic model, are valid here because the fatigue life is determined from the 
stresses in the elastic regime. However, fracture mechanisms are also of great interest in 
this area of study, in which case a plastic model should be developed that can accurately 
incorporate the elastic regime predicted here, as well as plastic behavior near the notch 
tip. Upon entering the plastic regime, however, future work should be aware of the many 
added factors that come into play, including strain hardening, crystal lattice rotation, 
creep, etc.. 

The orientations presented here have been studied by others, with varying degrees 
of success, yet few have conducted tests deviating from the <100> or <1 10> orientations; 
although, there is some degree of work done for the <1 1 1 > orientation. Since the elastic 
model’s simplicity makes it efficient, future research should entail testing for other 
orientations with a similar model, with the ultimate goal eventually being an accurate 
elastic-plastic model. Also, due to the variation in past research where sectors are 
determined by stress or strain fields, a study should be done to compare the stress and 
strain fields and the resulting sectors for the same load and model to confirm their 


correlation. 



APPENDIX A 

EXAMPLE COORDINATE TRANSFORMATION AND ACCURACY CHECKS 
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Checks for Accuracy 


All should equal zero: 
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All checks affirm a proper transformation. 
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dislocation) mechanisms, and no crystal lattice rotation. Finite element analysis is 
capable of accounting for changing temperature effects. However, in order to simplify 
the model and collaborate more closely with the MSE department, the FEA applied 
material properties at a constant (room) temperature. The exclusion of microstructure 
has also been discussed. The last assumption, however, is valid for this elastic model; 
crystal lattice rotation generally requires a few percent strain in a single slip system to 
occur (Stouffer and Dame, 1996) and the load used here is extremely low. The FEA here 
assumes only elastic deformation, and seeks only elastic stress and strain trends rather 
than specific values at fracture. Future research can easily build from the current model 
to incorporate plasticity, as well as creep or other desired behaviors. 



